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Abstract 

Recent theoretical investigations have unveiled a rich structure in the quantum chromodynamics 
(QCD) phase diagram which consists of quark gluon plasma (QGP) and the hadronic phases but 
also supports the existence of a cross-over transition ending at a critical end point (CEP). We 
find a too large variation in determination of the coordinates of the CEP in the temperature (T), 
baryon chemical potential (/x^) plane and, therefore, its identification in the current heavy-ion 
experiments becomes debatable. Here we use an equation of state (EOS) for a deconfined QGP 
using a thermodynamically consistent quasiparticle model involving quarks and gluons having 
thermal masses. We further use a thermodynamically consistent excluded volume model for the 
hadron gas (HG) which was recently proposed by us. Using these equations of state, a first order 
deconfining phase transition is constructed using Gibbs' criteria. This leads to an interesting 
finding that the phase transition line ends at a critical point (CEP) beyond which a cross-over 
region exists. Using our thermal HG model, we obtain a chemical freeze out curve and we find 
that the CEP lies in close proximity to this curve as proposed by some authors. The coordinates 
of CEP are found to lie within the reach of RHIC experiment. 
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I. INTRODUCTION 



QCD, the non-abelian gauge theory of strong interaction, predicts a phase transition be- 
tween a hot and dense hadron gas (HG) and the quark-gluon dominated phase which is 
called as quark gluon plasma (QGP) [1,2]. However, even after an intensive experimental as 
well as theoretical research spreading over the last two decades, our knowledge regarding the 
properties and signals of QGP is still very limited [3,4]. Even the phase boundary between 
the two phases remains in the literature as a conjectured one because the non-perturbative 
aspects of QCD are still dominant near the region close to the phase transition. In order 
to test the appearance of a deconfined QGP which exists in a transient phase, we need a 
proper understanding of its subsequent hadronization process about which our knowledge 
is really very poor. The numerical method of lattice QCD can properly describe both the 
phases i.e., QGP and HG. However, lattice QCD studies have yielded results for finite and 
large values of temperature T and /ib—O and now we have surmounted difficulties in getting 
results for small, non vanishing values oi /ib [5,6]. Thus we feel an urgent need to formulate a 
phenomenological model which can successfully reproduce the lattice QCD data and hence 
we can further use it to obtain the properties of QGP and to determine the QCD phase 
diagram in the entire (T, /x^) plane. 

In this paper, we present a thermodynamically self- consistent quasiparticle model of QCD 
which describes a gas of quasiparticles with effective masses generated through the inter- 
actions among its basic constituents. This model has been found to work well above and 
around critical temperature Tc- In order to describe the low energy HG phase, we work 
with our own model which has been found to give a proper description of the hot and dense 
HG [7]. The investigation of the structure of the QCD phase diagram is one of the most 
important and challenging topics in the nuclear and particle physics today. The precise 
determination of the phase boundary between QGP and HG at high temperature T and 
small /iB has been a subject of intense research in recent years from experimental as well as 
theoretical points of view [8-9] . Lattice simulations first revealed that the transition between 
HG and QGP phase at /jLb—O and large T is a crossover transition and there were further 
indications that the crossover transition turned into a first order chiral phase transition for 
nonvanishing and finite values of fiB [10]. Several attempts have since been made to locate 
precisely the critical end point (CEP) i.e., an ending point of the first order chiral transition 



as fj,B decreases [11]. Although the existence of CEP was predicted a long time ago by a 
few lattice calculations, the absence of the CEP in the phase diagram was also noticed in 
some recent lattice calculations [12,13]. Thus the location and the existence of the CEP 
in the phase diagram is still a matter of debate. Therefore, it is worthwhile to investigate 
the precise location of the CEP and to determine its properties in detail with the help of 
various phenomenological models [9]. However, confusion prevails since the coordinates of 
the CEP in the (T, hb) plane vary wildly in various models. Moreover, we are still not 
certain whether the conjectured phase boundary represents the chiral and/or deconfining 
phase transition line. Various calculations based on lattice QCD and/or effective models 
work with the basic assumption that the finite fis chiral phase transition is first order and 
hence the ending point of the line should give the CEP which is a second-order phase tran- 
sition point. However, we know that the chiral symmetry is broken in the colour flavour 
locked (CFL) region which has extremely large values ol jiB- Theoretically, this topic has 
largely been investigated using several phenomenological models which can give results in 
the entire (T, hb) plane, whereas ab initio calculations are still limited to very small values 
in ^B- However, all these efforts result in a very wide variation in the coordinates of the 
CEP [11]. In a recent paper, we proposed a new EOS for HG fireball where the geometrical 
size of the baryons in HG is explicitly incorporated as the excluded volume correction in a 
thermodynamically consistent manner [7]. Furthermore, we used a bag model EOS for the 
QGP phase and a first order phase transition is constructed by equating pressures of both 
the phases using Gibbs' criteria. We thus obtain an interesting result that such a simple 
picture not only reproduces the entire conjectured phase boundary but the first order phase 
transition line ends at a CEP beyond which a crossover region persists [7]. The cordinates 
of CEP as obtained in our calculation are found to be compatible with the prediction of a 
recent lattice gauge calculation of Gavai and Gupta [14]. Most importantly, we find here 
a deconfining phase transition in contrast to other calculations where the phase boundary 
depicts a chiral symmetry restoring phase transition. However, bag model is often a crude 
description to a gas of weakly interacting QGP and the nonperturbative effect in this model 
arises from the pressure of the vaccum through the use of a phenomenological bag constant. 
In the present work, we consider QGP as a system of quasiparticles which are quarks and 
gluons possessing temperature-dependent masses arising due to vacuum interactions [15-18]. 
Recently these models are made thermodynamically self consistent by incorporating suit- 



able corrections in two different ways and fience are referred as QPM I and QPM II in the 
following. Moreover, these descriptions were initially given for the gluon plasma only. We 
extend both the models for a QGP with finite baryon chemical potential /ib and we further 
assume their validity starting from a threshold temperature Tq = lOOMeV below the crit- 
ical temperature Tc and we adjust the parameters of the models accordingly. We further 
compare their predictions for the energy density e and the pressure p of the QGP with those 
obtained previously from the lattice simulation approach. Furthermore, we obtain a new 
EOS for the hot and dense HG as formulated in the previous paper and construct the phase 
boundary by equating the QGP pressure with that of HG pressure.We thus determine the 
critical parameters in the Gibbs' construction of a first order deconfining phase transition. 
We draw the phase boundary line and the end of the line determines the coordinates of 
the GEP. We also find the existence of a crossover transition lying beyond GEP. Finally we 
compare our findings with those obtained in various other models. 

The plan of the paper runs as follows. There are two types of quasiparticle models which 
are thermodynamically self consistent. In section II, we describe quasiparticles (QP), their 
corresponding equation of state and discuss the criterion of thermodynamical consistency. 
In section III, we describe the first thermodynamically consistent quasiparticle model [15] 
of Gorenstein and Yang, the dependence of the quark and gluon masses on the temperature 
and the extension of the model to describe the physics at finite /ib- In section IV, we discuss 
the second thermodynamically consistent quasiparticle model [16] of Bannur and its exten- 
sion for the finite /ib- In section V, we give our formulation for an excluded volume model 
for the hot and dense hadron gas, and we discuss about its thermodynamical consistency. 
Section VI summarizes our results and discussions. 

II. QUASIPARTICLE MODELS (QPM) 

Quasiparticles are the thermal excitations of the interacting quarks and gluons.The quasi- 
particle model in QGD is a phenomenological model which is widely used to describe the 
non-ideal bahaviour of QGP near the phase transition points. It was first proposed by 
Golviznin and Satz [17] and then by Peshier et. al. [18] to explain the EOS of QGP 
obtained from lattice gauge simulation of QGD at finite temperature. In quasiparticle mod- 
els, the system of interacting massless quarks and gluons can be effectively described as 



an ideal gas of 'massive' noninteracting quasiparticles. The mass of these quasiparticles is 
temperature-dependent and arises because of the interactions of quarks and gluons with the 
surrounding matter in the medium. These quasiparticles retain the quantum numbers of 

the real particles i.e., the quarks and gluons. It was assumed that energy uj and momentum 
k of the quasiparticles obey a simple dispersion relation : 

a;2(A;,T) = A;2 + m2(T), (1) 

where m(T) is the temperature-dependent mass of the quasiparticle. The pressure and 
energy density of the ideal gas of quasiparticles are dependent on uj and m(T) and are given 
by [15]: 
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where d represents the degeneracy factor for quarks and/or gluons. However, Gorenstein 
and Yang pointed out that this model involves a thermodynamical inconsistency [15]. It 
did not satisfy the fundamental thermodynamic relation: e{T) — T '^'^j^^ — p{T). So they 
reformulated the statistical mechanics for a system whose constituents follow a medium- 
dependent dispersion relation and the above inconsistency problem was handled by them 
by introducing a temperature-dependent vacuum energy term which effectively cancelled 
the inconsistent term. Alternatively, Bannur also pointed out the reason for the above 
thermodynamical inconsistency [16]. If the particle mass in the system is not constant 
and it depends upon the medium, then the relation used between the pressure and grand 
canonical partition function does not hold good. So one can start from the definitions of 
the energy density and the average particle number density in the grand canonical ensemble 
formalism and in this way a different thermodynamically consistent quasiparticle model for 
QGP can be obtained. 



III. FIRST QUASIPARTICLE MODEL (QPM I) 

Gorenstein and Yang initially formulated a thermodynamically-consistent quasiparticle 
model for a gluon plasma at /ib—O and later they extended it for the QGP having a finite 
value oi fiB- In this model, the effective mass of the gluon changes with T and //s as follows 



[15]: 

ml{T) = ^g\T)T^(^l + ^\, (4) 

where Nc represents the number of colours. We have also taken A^c = 3 in our calculation. 
And, 

/ 

Here Nf is the number of flavours of quarks and A*/ is the quark chemical potential belonging 
to the flavour f. Similarly the effective mass of a quark of flavour f changes with T and /ib 
as [19]: 

Here g^{T) is the QCD running coupling constant. We have taken the following form for 
g\T) [20-21]: 
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(7) 

Where Ay is the QCD scale-fixing parameter which characterizes the strength of the 
interaction. We have taken A^ = 115MeV in our calculation. Here parameter a is equal to 
^ [22]. 

After reformulating the statistical mechanics and incorporating the additional medium 
contribution, the pressure p and energy density e for a system of quasiparticles can be written 
in a thermodynamically consistent manner as follows [15,18]: 

p{T,m)^Pid-B{T,i,B): (8) 

e{T,m) = eid + B{T,iiB). (9) 

The first term on the right hand side of both the equations are the standard ideal gas 
expressions given by Eq (1) and (2), respectively. The second term represents the medium 
contribution: 

B{T,i,B)=lim ^, (10) 



where Eq is the vacuum energy in the absence of quasiparticle excitations or zero point 
energy. However, this energy term is not a constant but depends upon /ib and T. The 
B{T, hb) can be derived as follows [23]: 

^ „ d r ,^ r k^dk 1 , , 

B{T,i,b)^B,-— dT —-- 11 

47r^ JTa Jo ^ [exp {-^)\ T 1 

where Bq is the integration constant i.e., the value of -B(T, /ib) sitT — Tq. The expression for 
B{T, hb) in Eq. (11) and the forms of Eq. (8-10) give a compelling evidence that B{T, hb) 
may be treated as T and /is-dependent bag constant term [24]. In our calculation we take 
bI'^ = 185Me1/ and Tq = lOOMel/. We have also taken //^ = 3//^. 



IV. SECOND QUASIPARTICLE MODEL (QPM II) 



The second thermodynamically consistent quasiparticle model is given by Bannur [16]. 
Bannur first figured out why there exists a thermodynamical incosistency in the quasiparticle 
description of Peshier et. al. The relation between pressure and grand canonical partition 
function can not hold good if the particles of the system have a temperature-dependent mass. 
So he used the definition of average energy and average number of particles and derived all 
the thermodynamical quantities from them in a consistent manner. In this model, the 
effective mass of the gluon is the same as given in Eq (4). However, the effective mass of 
the quarks involves the following relation: 



T^l = "^go + V2mgomth + m^^. 



(12) 



were mgo is the rest mass of the quarks. In this calculation, we have used m^o = 8MeV 
for two light quarks {u,d), and rugo — SOMeV for strange quark. In the above Eq (12) nith 
represents the thermal mass of the quarks and it can be written as 



8Nr 



T2 + ^ 



(13) 



Taking these values for the effective masses, energy density can be derived from the grand 
canonical partition function in a thermodynamically consistent manner and is given as : 



OO ^ r , , ■ 

E /I + i-ly-%cosh{^^/T)eixgl) + (-l)'-i^e,(xj) 



(14) 



with ei{xil) — {xil)^Ki{xil) + 3{xiiyK2{xil), where Ki and K2 are the modified Bessel 
functions with Xi — and index i runs for gluons, up-down quarks q, and strange quark s. 
Here di are the degeneracies associated with the internal degrees of freedom. Now, by using 
the thermodynamic relation e = — p, pressure of system at /^^ = can be obtained as: 

T To Jt, r2 ' y J 

where po is the pressure at a reference temperature Tq. We have used po=0 at To=100 
MeV in our calculation. Using the relation between the number density Ug and the grand 
canonical partition function, we can get the pressure for a system at finite hb '■ 

P{T,Hq) = P{T,0) + Uqd^q. (16) 

Jo 

Thus all the thermodynamical quantities can be obtained in a consistent way by using this 
model. 



V. EOS FOR A HADRON GAS 



There is no deconfinement transition, if the hadron gas consists of point-like particles, 
and consequently HG pressure is always larger than QGP pressure. Therefore, inclusion of a 
repulsive interaction between two baryons having a hard-core size reduces the HG pressure 
and hence it stabilizes the formation of QGP at high baryon densities. Recently we have 
proposed a thcrmodynamically consistent excluded volume model for hot and dense hadron 
gas (HG). In this model, the grand canonical partition function for the HG with full quantum 
statistics and after incorporating excluded volume correction can be written as [25] : 
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where gi is the degeneracy factor of ith species of baryons, is the energy of the particle 
{Ei — ■\/k'^ -\- ml), is the eigenvolume of one baryon of ith species and Ylij ^j^f 
total occupied volume and Nj represents total number of baryons of jth species. 
Now we can write Eq.(17) as: 



InZr ^V{l-J2 nTVf)hK (18) 

3 

where /j represents the integral: 

i.^^rj^ I (19) 

and Aj = exp(^) is the fugacity of the particle, is the number density of jth type of 
baryons after excluded volume correction and can be obtained from Eq.(18) as: 



This leads to a transcendental equation as 



-'I- R)IX - liX-g^ + A,'(l - R)l'i (21) 

where /■ is the partial derivative of /j with respect to Ajand R — J2i '^T^i fractional 
occupied volume. We can write R in an operator equation as follows [7]: 

R^Ri + (lR (22) 

where Ri — with i?° = ^ T^V^i + ^'i^i^h ''^i density of pointlike baryons of 

ith species and the operator Cl has the form : 

" = -T^E"?V°A.|- (23) 

i 

Using Neumann iteration method and retaining the series upto Cl^ term, we get 

R^ Ri + nRi + n^Ri (24) 

Eq.(24) can be solved numerically. Finally, we get for the total pressure [25] of the hadron 
gas: 

i i 

In (25), the first term represents the pressure due to all types of baryons where excluded 
volume correction is incorporated and the second term gives the total pressure due to all 
mesons in HG having a pointlike size. This makes it clear that we consider the hard-core 
repulsion arising between two baryons which possess a hard-core size. In this calculation. 



dR 



we have taken an equal volume — for each type of baryon with a hard-core radius 
r = 0.8 fm. We have taken all baryons and mesons and their resonances having masses upto 
2GeV/c^ in our calculation for HG pressure. Wc have also used the condition of strangeness 
conservation by putting Si{nf — nf) = 0, where Si is the strangeness quantum number of 
the ith hadron, and is the strange (anti-strange) hadron density, respectively. Using 

this constraint equation, we get the value of strange chemical potential in terms oi /ib- We 
want to stress here that the form of this model used under Boltzmann approximation has 
been found to describe [26-27] very well the observed multiplicities and the ratios of the 
particles in heavy-ion collisions. 



VI. RESULTS AND DISCUSSION 



In order to demonstrate that both types of quasiparticle models reproduce the lattice 
results with the value of parameters chosen here, we show in Fig. 1, the results of our 
calculations for the variations of energy density with respect to temperature at different /ib- 
We find that the predictions from both these models (QPM I and QPM II) compare well with 
the lattice data [28] , although it cannot be regarded as a very good fit to the data. However, 
the results of QPM may improve if wc use thermal gluon mass as m^(T) = gf^(T)T^/3 
instead of rri^lT) — g'^{T)T'^/2 as pointed out by Bannur [16,21]. 

In Fig. 2, we have presented the results of our calculations for the QGP pressure p/T^ in 
both the quasiparticle models and shown its temperature variation at different values oi /ib- 
We compare our results with those recently reported in lattice simulations [28,29]. We find 
that the fits by QPM II look much better than those given by QPM I. These results give us 
extreme confidence in both types of quasiparticle models being used as phenomenological 
models. Although phenomenology cannot work as a substitute for a formal theory like QCD. 
However, since the utihty of lattice QCD calculations at very large hb is still not possible, 
we take the help of quasiparticle model which has a few adjustable parameters. It is now 
widely used to describe the non-ideal behaviour of QGP observed near the critical line. 
Moreover, we attempt to extend its uses at lower values of temperature, eg., Tq < (we take 
To = lOOMeV). Usually authors have studied the quasiparticle models above Tc only and 
therefore, the rapid rise of pressure and energy density at or around T^. is not properly taken 
care of in these models. Our method of obtaining critical parameters (Tc, /ic) involves the use 
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FIG. 1: The calculation of energy density of QGP in QPM I (by a dashed curve at /xb/T=0.0 
and by a dotted curve at //b/T=2.0) and in QPM II (by a sohd hne at //b/T=0.0 and by a 
dashed-dotted curve at //s/T=2.0). Square points are lattice data at /is/T=0.0 [28]. 

of Gibbs' equilibrium criteria of equating HG and QGP pressures and to determine where 
these pressure lines intersect each other. Therefore, we want to know precisely the values 
of QGP pressure at temperatures T'o(< Tc) also. The comparison of our calculations with 
the lattice results yields the required test about the suitability of quasiparticle EOS for QGP. 

In Fig.3, we have shown the phase boundary obtained in our model. Surprisingly we 
again find here that the first-order deconfining phase transition line ends at a critical end 
point (CEP) and the coordinates of CEP are (T^ = 183MeV, jj^ = IQQMeV) in QPM I and 
{Tc — 166MeV, //c = IbbMeV) in QPM II. It is interesting to find that the critical points 
obtained by us he closer to CEP of some lattice calculation [14] . These points also compare 
well with the coordinates of CEP (T^ = 160MeV, /jLc = 156MeV) obtained by us in the 
previous publication [7] using bag model calculation. We also find a crossover region existing 
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FIG. 2: The variation of QGP pressure with the temperature at different ijlb- In QPM I we have 
shown the results by dashed curve at //b/T=0.0 and by a dotted curve at //b/T=2.0. Similarly in 
QPM II, the solid line represents results at ;^b/T=0.0 and the dashed-dotted curve at /^b/T=2.0, 
respectively. Square and triangular points are lattice data at /xb/T=0.0 and at /xb/T=2.0, respec- 
tively [28,29]. 

beyond the critical point where HG pressure which is solely dominated by mesonic pressure 
term in Eq.(25), is always less than the QGP pressure. Therefore, no phase transition exists 
in this region. Since the temperature is much higher, the thermal fluctuations break mesonic 
constituents of HG into quarks, antiquarks and gluons. We have tabulated in Table I the 
location of CEP obtained from various calculations for a comparison. For convenience we 
have shown above the dark solid line, all the values obtained in SU(3) flavour calculations. 
Below the solid hue, the values of SU(2) are also shown in order to make the comparison 
complete. 

We find that there exists a very wide variation in the coordinates of GEP obtained in 
different models. We notice that critical end points obtained in deconfining phase transition 




FIG. 3: The location of QCD critical point in QCD phase diagram. Pi is the phase boundary in 
bag model (BM), P2 is the phase boundary in QPM I and P3 is the phase boundary in QPM II. 
Fi is the chemical freezeout line obtained using our HG model. Ci{Tc = ISSMeV, /ic = IQQMeV) 
is the CEP on P2 obtained in QPM I and C2(Tc = 166Mey,/ic = ISSMeV) is the CEP on P3 
obtained in QPM II. The labels used in the figure are explained in the table I. 



are usually located at jiB < 200MeV whereas chiral CEP have much larger /x^. So there 
is a good chance for observing CEP at RHIC by using energy scan [47]. We also find that 
CEP obtained in our models almost overlaps with the points on the freezeout curve. The 
freezeout point occurs close to RHIC energy and hence the fluctuations in multiplicity etc. 
can experimentally provide a clear signal for CEP. 

Before we conclude we would like to discuss the recent lattice data obtained by de For- 
crand and Philipsen [12, 13] who used 2+1 and 3 flavors staggered fermions and a taylor 
expansion in fiq/T to study the curvature of the critical surface at very light quark masses 
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TABLE I: Coordinates of CEP obtained in different models. The last column gives the correspond- 
ing label used in Fig. 3 and Reference [11] 

close to /ig = surface. They noticed that the critical surface bends so that the first order 
region shrinks at higher quark masses and hence they conclude that there is no critical point 
at finite chemical potential. However, it is speculated that the critical surface bends back at 
larger fig and the critical point may again reappear. In fact a recent NJL model calculation 
lends support to this speculation [48]. 

In summary we have demonstrated the occurrence of CEP in a deconfining first order 
phase transition constructed by using quasiparticle model for QGP and a new thermody- 



namically consistent excluded volume model for hot and dense HG. In our model, we assign 
a hard core volume for baryons and mesons are treated as pointlike particles. So at higher 
temperatures, mesons can fuse into one another but baryons occupy space. As increases, 
we find that the fractional occupied volume R by baryons increases and hence the mobility 
of baryons decreases fast. The physical mechanism in our model is similar to the percolation 
model [49] where a first order phase transition results due to "jamming" of baryons in the 
HG. Thus our finding lends support to the idea of realising a phase transition by modelling 
the interactions existing in the HG in a suitable way. 
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